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ABSTRACT 

The far-field radiation of turbofan noise from the inlet of an engine 
has recently received considerable study. Relatively simple approximate ex- 
pressions have been developed for single mode radiation as a function of cut- 
off ratio and for multimodal radiation. Convective effects upon radiation 
which are essential for interpreting static engine tests have been determined. 
This information is essential for the design of Inlet suppressors. This paper 
represents an initial effort to provide similar approximate expressions for 
the noise radiation from the aft duct. The results of approximate aft radia- 
tion equations compare favorably to more exact Wiener-Hopf radiation results. 
Refraction as well as convective effects in the multiple flow streams must be 
considered. The peak in the radiation pattern, which occurs nearly at engine 
sideline, is composed of modes with relatively large cut-off ratios. This im- 
plies that aft fan radiation will be inherently more difficult to suppress 
than the fan inlet noise. The theoretical multimodal radiation pattern is 
compared to experimental data for the first two harmonics of blade passage 
frequency for three full scale fans at two speeds. The agreement between 
theory and experiment is quite good. 


INTRODUCTION 


Far-field radiation theory for sound radiation from engine ducts is essen- 
tial to noise suppressor studies. First, it is required to translate the 
sound attenuation occurring in the duct to the more useful information of 
sound attenuation in the far-field 1 *^ of the engine. Second, the radiation 
theory, when compared to noise directivity data, can be used to infer informa- 
tion about either the modal power distributions*^*^ or the acoustic power dis- 
tribution versus mode cut-off ratio® which exists in the engine duct. 

Sound radiation from the inlet duct has been studied extensively. Tyler 
and Sofrin® developed a flanged duct radiation model which was shown to pro- 
vide a reasonable representation to experimental data. Lansing, et al.9 and 
Homicz and Lordi^® developed Wiener-Hopf radiation solutions for an unflanged 
duct. Choli presented sound radiation solutions for hyperboloidal inlet ducts 
which were favorably compared to experimental data by Clark and Cho. 12 There 
are serious limitations on all of these exact analytical solutions as presented, 
The flanged duct and hyperboloidal inlet solutions are currently formulated 
only for zero steady flow both in the duct ana in the surrounding environment. 
The Wiener-Hopf solutions can accomodate steady flow but only for the same 
flow velocity in the duct and in the surrounding region, or else for a tube of 
high velocity flow separated from the ambient flow by a slip layer. Neither 
of these flow situations are realistic for real engine inlets, especially at 
static test conditions. Under static test, the flow velocity is very high in 
the duct while outside there is essentially zero flow at all angles from the 
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inlet axis. This more realistic but complicated flow situation was analyzed 
approximately by Rice, et al.4 For the two cases, no flow and equal flow 
everywhere, which can be handled by the Wiener-Hopf method, the far-field 
radiation angle for the peak of the principal lobe of radiation was found to 
be coincident with the angle of propagation with respect to the axis for the 
mode in the inlet duct. This observation led to using simple angle trans- 
formations to account for the local effect of flow on the principal lobe of 
radiation in going from the high Mach number in the duct to a different flow 
Mach number in the ambient medium. This method was indirectly confirmed by 
its ability to explain far-field directivity attenuation patterns 1 * 2 and more 
directly, by observation of single mode direcitivlty patterns in a static en- 
gine test. 5 littidmann, et al. 5 also generalized the above approximate theory 
by extending the duct flow transformation to include any angle instead of 
only the angle of the principal lobe of radiation considered in reference 4. 

Some other observations concerning the far-field radiation theories 
should be made. Lansing, et r 1.9 have shown that except in the vicinity of 
the duct flange, the flanged and unflanged duct radiation solutions are almost 
identical, Homicz and LordiiO have amplified this correspondence by deriving 
the zeroes (cancellations) in the unflanged duct Wiener-Hopf solutions and 
they are seen to coincide with those of the flanged duct solutions. Flanged 
duct radiation solutions are much simpler in form and use than the Wiener-Hopf 
solutions, which require complex computer programs for evaluation. The sim- 
pler closed form solutions should always be used when they are sufficiently 
accurate. They also yield to approximations which allow the direct derivation 
of simple expressions for multimodal radiation patterns. These multimodal 
radiation patterns have been presented by Rice^ for zero steady flow and have 
been found to adequately reproduce exact radiation calculations . 6 These 
simple expressions have been found useful in inferring the distribution of 
acoustic power versus mode cut-off ratio for sound with randomly generated 
modes, such as broad-band noise or even tones found in static engine tests. ® 

The far-field radiation theory for the aft duct case has received much 
less attention than that for the inlet. This may be due to the complex multi- 
ple flow streams at several Mach numbers which make up a typical exhaust duct 
system. However, compared to the inlet case, the Wiener-Hopf solution should 
provide a more realistic modelling of sound propagation through the exhaust 
flow streams with slip layers. Wiener-Hopf radiation solutions have been pre- 
sented by Savkar and Edelfelt 15 and Savkar 1 '* for a single jet stream with a 
flow mismatch with the ambient environment. These exact solutions are used 
as a base for checking simpler approximate solutions. 

The purpose for the work reported in this paper is thus to generate ap- 
proximate radiation expressions for the aft duct, in the same spirit as was 
done for the inlet, and check these simple expressions against the available 
exact solutions. The starting point for this approximate analysis is the 
known propagation angle for a mode in the exhaust duct with flow. An angle 
transformation is performed to obtain the direction of the vector normal to 
the wavefront. Refraction is then considered through the slip layer and, 
finally, another angle transformation is required if there is an ambient flow 
velocity. Also, the flanged duct radiation equation is cast in a form which 
is usable for the aft duct. The above transformations are required to generate 
the far-field radiation pattern. The simplified radiation patterns are then 
compared to the exact single mode solutions. A simple multimodal radiation 
pattern is then generated from the approximate expressions. Results from this 
simple equation and data from three full scale fans at two frequencies and 
speeds are compared to a laboriously obtained exact multimodal radiation pattern. 
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The approximate expressions obtained here are essential for the develop- 
ment of an aft duct liner design and evaluation program analogous to the inlet 
program reported in reference 1. 


SYMBOLS 

c speed of sound, m/sec 

c„ c in exhaust duct and jet, m/sec 

c^ c in fiir-field surrounding medium, m/sec 

D duct diameter, m 

f frequency, Hz 

i far- field angle index in equation (38) 

j' derivative of Bessel function of first kind of order m with respect 

m to its argument 

Mp Mach number in exhaust duct and jet 

M, Mach number in surrounding medium 

t 

m spinning mode lobe number 

N relative number of modes, see equations (34) and (35) 

n . 

P acoustic pressure, N/m* 

fch 2 

P, acoustic pressure at i c far- field angle, M/m 

* 2 
P £w multimodal P in far-field with duct flow, see equation (36), N/m 

tM 

P. principal lobe peak P in far-field with duct flow, see equation 

£P (31), N/m2 

P_ principal lobe peak P, flanged duct, no flow, see equation (22), 

P N/m2 

R radius from duct exit to far-field observation point, m 

£S n far-field area subtended by principal lobe, no flow, see equations 

0 (26) and (29), n.2 

far-field area subtended by principal lobe, with duct flow, see equa- 
tion (30), m2 

SPL^ multimodal sound pressure level in far-field with duct flow 

W acoustic power, Nm/sec 

a hardwall eigenvalue for (m, (-0 mode 

m|a _____ 

p f/T- Mq 

q frequency parameter, fD/c 

q frequency parameter transformed to steady flow, q/p 

|i,v radial mode index, p = 1 is lowest order 

Ip mode cut-off ratio with duct flow, see equation (8) 

see equation (32) 
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fz 

fzs 


tfr 


OP 


ifr 


ft 
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02 


see equation (33) 

mode cut-off ratio, no flow 
3 

density, kg/m 

angle between phase velocity vector and duct axis, deg 
9 in duct with flow, deg 
9 in far- field, deg 

angle between group velocity vector and duct axis, deg 

group velocity angle transformation to introduce steady flow, see 
equation (5), deg 

ft in duct with flow, deg 

$ in far-field with duct flow, deg 

ft^ at principal lobe peak, deg 

tg at zeroes of radiation pattern, deg 

fg at boundary of zone of silence, deg 

t in far-field, no flow, deg 

at principal lobe peak, deg 

$ at zero of principal lobe below peak, deg 

ft at zero of principal lobe above peak, deg 


APPROXIMATE RADIATION THEORY 

In this section the flanged duct radiation equations which are to be used 
are presented. The several angle transformations required are also given. The 
geometry considered is a circular jet with a flow Mach number M D discharging 
into a surrounding medium with Mach number Mf (symbols are defined in the 
Symbol List). The subscript D refers to the duct and ItB related jet while 
subscript f refers to the far-field or surrounding medium. The angles 9 
are between the duct axis and the phase velocity vector (normal to the wave 
front) while angles ft are between the duct axis and the group velocity vec- 
tor. For zero flow 9 and ft are the same. Although the transformation are 
derived for the general condition, the surrounding medium will be assumed at 
rest as in a static engine test for the results discussed in this paper. 


Radiation Equations and Angle Transformations 


The far-field mean-squared pressure radiation expression for a flanged 
circular duct for equal acoustic power per mode3»® and zero flow can be ex- 
pressed as, 


2 V" 1 '(^f) ■ i ” 2 *0 6in < '0> ]2 

P 3 — 1 2 


/ w 9 


r 2 i 

i-(f-) 


(•*) - 

\m/ 


\ 7U| / V 


( 1 ) 
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whore the subscript on % denotes zero velocity and the frequency parameter 
is 


n 


fD 

0 - 


( 2 ) 


and the mode eigenvalue is determine from 

j ;<v> ■ 0 < 3 > 

An approximate version of equation (1) has been reported^ using the sine ap- 
proximation for the Bessel function and requiring (m/omy)* « 1. This expres 
aion, which will be useful in a later discussion, is, 



The first transformation required is to introduce uniform flow into equa- 
tion (1) by using a coordinate stretching technique. 9,10,15 mfl y seem physi- 
cally strange to have uniform steady flow with a flanged duct configuration, 
but the final results will justify its use. The transformation involves re- 
placing ¥q by ¥ such that 


sin ¥ « 


p sin 

|/ 1 - sin 2 f D 


and replacing q by 



where 


p = y'l - m£ 

In addition the mode cut-off ratio with flow is, 



(5) 


( 6 ) 


(7) 


( 3 ) 


When equations (5) to (8) are used in equation (1) the result is, 
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Equation (9) describes the radiation in a uniform flow field and should ade- 
quately reproduce the results for an uttflanged duct, except in the vicinity 
of the flange where some error will occv--. With the Mach number M D in the 
radiation expression, there are more propagating modes than with Mp * 0. 

This can be seen from equation (8), where for a given frequency parameter (q) 
and duct mode (eigenvalue a^) the cut-off ratio is higher with uniform flow. 
Cut-off occurs when Ip * 1. 

The second transformation yields the phase velocity vector propagation 
angle q>p, which is required for the refraction expression through the flow 
discontinuity at the jet boundary. Figure 1 shows the several angles involved. 
The angle tp D is related to the group velocity vector angle fp by, 


COS 


Mp + cos cp D 


+ m£ + 2^ cos cp D 


( 10 ) 


as discussed in reference 4. Equation (10) is then Inverted to yield, 


cos 


(p D * -Mjj sin 2 f D + sin 4 f D - (l + M^sin 2 ^ + 1 




>in 2 t D + cos f D |/ 1 - M^ sin 2 t D 


(U) 


The refraction through the slip layer between the aft duct jet and the 
surrounding medium can now be calculated. This is done in the simplest way by 
considering the refraction of an oblique plane wave upon a plane slip layer, 
as can be found in reference 16, The propagation angle of the phase velocity 
emerging from the slip layer into the surrounding medium is, 


cos cp f 


C f COS 

C D + <V*D " c f M f> cos * D 


( 12 ) 


Finally, the direction of the group velocity vector in the far- field can be 
derived in a manner parallel to that of equation (10) to obtain, 

+ cos cp. 

COS 

f^l + M 2 + 2M f cos cp f 


( 13 ) 
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Thu angle Is the far- field radiation angle used for plotting the trans- 
formed flanged duct radiation pattern. 

The simpler case will be considered in all that follows, in which die 
surrounding medium is static (Mf * 0) and both media have the same temperature 
so that Cf * cp * c. Equations (12) and (13) then simplify to, 


cos <p 

00fi ■ 1 + M„ cos tp D 


and 


cos fg *• cos <p^ 


(14) 


(15) 


The angle transformations from equations (11), (14), and (15) can be combined 
to yield, 


-M d sin 2 t D + /m 2 sin 4 t D - (l + Mj)sin 2 * D + 1 

COS #£ * — — 

1 - sln 2 ^ D + Mp sinV D - (l + M^ain 2 ^ + 1 

-»B “l-^D + °° S *D » /l - '4 “ ln2 "’'D 

* (lb) 

1 - M^ sin 2 ^ + ^ cos f D - m£ sin 2 ^ 

Values of can now be inserted into equation (9) to calculate from 

which the sound pressure level can be calculated, and finally the radiation 
angle in the far-field can be calculated from equation (16). 


Comparison of Approximate and Exact Single Mode Radiation Theory 

The approximate radiation expressions as presented in the previous sec- 
tion were exercised for the (8, 4) mode (i.e., m ** 8, p « 4). This mode was 
chosen because it has several sidelobes on each side of the principal lobe. 

The convention used here is that the lowest radial mode has |i = 1. The exact 
Wiener-Hopf solution*^ was generated using a program available in the NASA- 
Langley computer. Both solutions are plotted in figure 2. The agreement be- 
tween the two curves is very good. Most important is the fact that the prin- 
cipal lobe location and shape are very accurately reproduced. The zeroes or 
nodes are also accurately located. Some error occurs in the amplitude of the 
side lobes . The accuracy of the approximate method is all the more surprising 
since there are several other physical phenomena (included in the exact solu- 
tion) which could be added to the theory but have not been completed at this 

time. Internal reflection at the slip layer as a function of propagation 
angle should be included. The various angles involved in the flanged duct 
solution are shifted differently by the refraction at the slip layer. This 
involves far-field area changes which should be accounted for (some discussion 
occurs in a later section) to retain acoustic power equality. Also it is 

known that at the flange, the flanged duct solution yields results which are 
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high compared to the unflanged duct solution, due to reflection off the flange. 
This could be corrected approximately, 


Location of the Principal Lobe and the Zeroes 

The critical values of the approximate far-field radiation theory will 
occur when the argument of the derivative of the Bessel function in equation 
(9) is equal to a duct eigenvalue, i.e., 


nf) sin 


a 


mV 


(17) 


D 


If the (m, u) mode is considered, then the location of the peak of the prin- 
cipal lobe will be found when v « u. Note the denominator of equation (9) 
also vanishes at this point and a finite limit is found for the amplitude. 
When v ^ |4 only the numerator of equation (9) vanishes and the several 
zeroes are located. By using the transformations given by equations (11), 
(14), and (15), equation (17) can be reduced to, 


cos fgp sw 


-“d 


(«l\ J(jsl\ . 

vW/ f \ a mn/ 


mU w \ m|i, 

which gives the principal lobe peak location and, 






(18) 


cos * fz « 


““d 


(JB.) | 
vmv/ V \ a mv/ 


(t) 


(19) 


I— -J 

which gives the zeroes for the (m, u) mode. Note that the mode cut-off ratio 
(S D ) for the (m, u) mode can be introduced into equation (18) to yield, 



-Mp + P 2& D 

P ^d(^d + **0 ✓w) 


( 20 ) 


Note that the location of the principal lobe is a function only of the mode 
cut-off ratio and the Mach number. All modes with similar cut-off ratios will 
radiate similarly to the far-field just as in the zero Mach number case.^ 

The results of equations (18) and (19) are compared to exact Wiener-Hopf 
solutions in figure 3 for four modes. The approximate principal lobe peak and 
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che zeroes are denoted by arrows, Note that the agreement Is excellent, 

Some derivation has been noted In other plots (not shown) when the principal 
lobe is located far ahead of the aideline (#£ 90°), It is not known whether 

the approximate theory or the exact theory has deficiencies in this range. 

It should be noted that the equations in this section were derived for a 
static surrounding medium. For Mf ^ 0, equations (12) and (13) should be 
used in the derivation Instead of equations (14) and (15). 


Zone of Relative Silence 

If * 0 is used in equation (16) the limiting angle for the zone of 
silence can be defined as, 


cos 



1 + 


>*0 


( 21 ) 


For the example radiation patterns shown so far in this paper with Mp * 0.6, 
the boundary of the zone of silence is at ff * 51.32°* Thus, in the approxi- 
mate theory used here, there is no radiation at all predicted between the aft 
duct axis and tf » 51,32°. This result is of course the consequence of the 
simple refraction model (equations (12) or (14)) used across the slip layer. 
The exact Wiener-Hopf technique also shows an abrupt change at tf * 51.3° as 
shown in figures 2 and 3. However, these exact solutions show some radiation 
within the zone of silence although the levels are very low. 


Multimodal Far-Field Radiation Model 

In this section, the approximate multimodal radiation equation developed 
in reference 6 will be extended to include the duct jet flow and the refrac- 
tion through the slip layer. This simple approximate equation will then be 
compared to the exact multimodal Wiener-Hopf solution, which is quite cumber- 
some to generate. Equal acoustic power per mode will be considered here, but 
the development could be easily generalized by using a cut-off ratio biasing 
function. ^ 

Approximate theory . - The technique used here will be illustrated by ex- 
panding upon the zero Mach number theory developed in reference 6. The ap- 
proximate single mode radiation equation (4) will be used as a starting point 
since it is quite a good approximation to equation (1) but is much easier to 
use. At the peak of the principal lobe in the far-field, equation (4) yields, 6 



T) COS f. 


OP 


2 sin % 


OP 


( 22 ) 


where the zero Mach number cut-off ratio is obtained from equation (8) with 
Mj) = 0 and the principal lobe location is given by, 6 


sin f, 


OP 


_ 1 _ 

g 0 


(23) 
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From equation (4) it can be seen that the zeroes in P^ extend around the 
principal lobe peak at angle* 


sin 




1 

n 


sin 


A 

n 


(24) 


sin #Q 2 w 



sin ?Qp 



(25) 


The far-f,ield incremental area subtended by the principal lobe can thus be 
estimated by, 


£S 


0 



sin * Q df Q n 2nR 2 sin t 0p (t 02 



(26) 


where a small angle increment is assumed, A first o*uer expansion of foi 
and fp2 around top can ma de assuming 1/q is very small so that 
to2 " ^01 can be estimated as, 



to 




Using equation (27) in equation (26) yields 

4nR 2 sin t 0p 
^0 * q cos t Qp 


(27) 


(2B) 


Finally the acoustic power in the principal lobe can be estimated by, 



Thus, at a constant radius the acoustic power in the principal lobe is constant 
no matter where it is located. The far- field pressure fall-off toward the 
side-line is due only to an Increase in the incremental far- field area. The 
R z factor remains in equation (29) since it was eliminated in the derivation 
of the directivity factor presented by equations (1) or (4). 

In a manner very similar to the above, the far- field area subtended by 
the principal lobe of radiation from a duct with jet flow and a surrounding 
medium separated by a slip layer can be derived. The derivation is far more 
complicated since the several angular transformations going from fp to 
to <p£ to must be used. The final result is, 


4nR 2 
£S £ « 


/T^Mp coa » fp - (l - >4)cos 2 » fp 

"["d + ( l • M S)“ 8 *fl>] 


( 30 ) 


u 


By analogy with equation (29) the principal lobe peak can be estimated by, 



2nR‘ 

48 



(3D 


As in reference 6 the mutlimodal pressure squared will be estimated by multi- 
plying ;he principal lobe peak by the number of modes whose principal lobe 
con contribute at a particular angle, By using the transformations of equa- 
tions (5) and (6) in equation (4), the cut-off ratio range which can contri- 
bute is seen to be from, 

, 6 sin ft* 

.JL * p p . a 

|/ 1 - flln 2 ft D ^ 

up to 

_L., g . ,1 l!a._. + a 

5 “ 2 fx »£ sin 2 # D 11 

These can be used in the integration of the modal density function*-? to obtain, 

6 D1 9 

9 4jT sin ft_ 

4 dS D « (34) 

& D q V 1 m Bitl \ 

D2 



(32) 


(33) 


This can be transformed to (using equations (10), (14) and (15)), 



The multimodal (equal acoustic power per mode) radiation pattern is then esti- 
mated by, 

4 ' 4" - 2 ( l • $ [“d + ( l - ' 4 )°°° **] < 36 > 

where the general far- field angle ff replaces ftfp. Note that when Mp * 0 
equation (36) reduces to 2 cos ftf, which is the same result as in reference 6. 
Equation (36) can be used only outside of the zone of silence since both N 
and 4Sf -♦ 0 as cos ftf -* 1/(1 + Mp), and P^ can be seen to peak at this 
angle. If this value of p|„ is used as a normalization, a sound pressure 
level can be calculated as, 
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SPL f 


** 10 Log 




(37) 


which can be seen to be a very simple expression for the multimodal, equal 
acoustic power per mode directivity, this expression will be compared to 
exact results in the next section. 

Exact theory . - The more exact multimodal (equal acoustic power per mode) 
results were generated using calculations performed on a NASA Langley computer 
with a Wlener-Ilopf radiation program . 13 only single mode results were avail- 
able, examples of which have been shown in figures 2 and 3. The condition for 
the calculations were M D * 0.6 and q » 7.111. At these conditions, it was 
estimated there would be 100 propagating modes. 17 The actual count showed 
that 101 modes could propagate. This provided enough modes to give a fairly 
smooth radiation pattern while keeping the tedious modal superposition within 
manageable limits. The approximate equation (37) shows that the radiation 
pattern for equal acoustic power per mode is a function only of duct Mach num- 
ber. 

Unfortunately, the exact single mode calculations gave no acoustic power 
number such that they could be readily combined. All results were given with 
the peak in the radiation pattern having a sound pressure level of 100 decibels. 
All 101 radiation patterns had to be normalized manually to obtain equal 
acoustic power per mode. Due to these time consuming procedures only one 
multimodal example will be shown. 

For each single mode radiation pattern the pressure squared tern, 
(arbitrarily normalized) was given at each 5° Increment in angle ($£). These 
values were area weighted to obtain a value which was proportional to the 
arbitrary acoustic power of the particular mode, using the equation (18), 






(38) 


where « 5°, Each value of p| was then normalized by the total power 
using equation (38), thus guaranteeing a particular total acoustic power for 
the mode. This was repeated for each of the 101 modes. The newly normalized 
value of p| for each of these modes was then added over all modes for a given 
angle (fixed i). This summation of pressure squared terms was then converted 
into a relative sound pressure level at each angle. The results, representing 
the combined effects of 101 modes all having equal acoustic power, ore Bhown 
in figure 4 as the circle symbols. The curve is relatively smooth with the 
peak occurring just forward of the zone of relative silence. The fall-off 
toward the inlet axis (180°) is seen to be quite gradual. A large reduction 
in sound pressure level occurs in the zone of silence just as in the single 
mode cases. 

The approximate radiation pattern for equal acoustic power per mod's, equa- 
tion (37), is also shown in figure 4. The agreement is good over most angles. 
The exceptions are at the extreme angles, where approximate theory is slightly 
high at the edge of the zone of silence and quite low toward the inlet axis 
(large fg). As mentioned previously, the approximate theory could be improved 
but it is probably good enough for most liner design purposes. Also note that 
in the zone of silence the approximate theory predicts no radiation while the 
more exact results show some radiation, although at considerably reduced level 
as compared to the peak. 
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Comparison with experimental data , - Noise data for three full scale fans 
tested statically without inflow control screens will be compared to the multi- 
modal theoretical results, Results for the blade passage tone and the second 
harmonic will be presented for both take-off and approach speeds. Tables I 
and II, which were obtained from information in references 19 to 22, present 
the design and aeroacoustic parameters for the fans designated A, B, and C. 

The noise data represent narrow band data from which the broadband noise has 
been subtracted as in reference 23. The theoretical results are from the 
more exact calculations with equal acoustic power per mode for Mp * 0.6 as 
shown in figure 4. The approximate theory, equation (37), showed about the 
same calculated results for the Mach number variations encountered here. 

Thus, the single theoretical curve will suffice for the qualitative compari- 
sons to be made here, It should be noted also that if data compares well with 
an equal acoustic power per mode radiation theory then a random noise genera- 
tion process is implied such as might result from blade-turbulence interactions. 
This appeared to be the case for static fan inlet data as shown in reference 6. 
The noise data for the blade passage fundamental should be cut-off if caused 
by rotor stator interaction. The second harmonic has rotor stator interaction 
modes which can propagate. 

The experimental data for the blade passage frequency fundamental at ap- 
proach speed are compared to theory in figure 5 for the three fans. The data 
for fans A and C show several similarities with the theory. For angles forward 
(larger angles) of the zone of silence the slope of the data is quite flat as 
is the theory. The deviation in the data here may be either scatter or domi- 
nant mode peaks. The data also show an abrupt fall-off in the zone of silence. 
Fan B conversely shows a more gradual fall-off into the zone of silence. No 
explanation is readily available for this behavior. 

In figure 6 for the blade passage frequency at take-off fan speed the 
theory-data comparison is very similar to figure 5. The higher Mach number, 
see table II, would cause the zone of silence to move to higher angles. Some 
hint of this may be seen in the data for fan A. The theoretical calculation 
has not been adjusted for higher Mp. 

The second harmonic noise data at take-off speed are compared to the the- 
ory in figure 7. These data are seen to be considerably different from those 
for the fundamental. Considerable dominant mode lobe patterns are observed 
for fans A and B. These might be expected since propagating modes are avail- 
able from the rotor-stator interaction. The considerable penetration of these 
dominant modes into the zone of silence cannot be explained at this time. 


ATTENUATION POTENTIAL FOR INLET AND EXHAUST SUPPRESSORS 

During the course of this study of radiation from an exhaust duct, a dis- 
tinct difference between an inlet and an exhaust duct was observed regarding 
the relationship between angle of radiation and mode cut-off ratio. The follow- 
ing initial observations are offered although the complete ramifications of the 
observations have not been investigated. Some additional factors which must be 
considered are discussed at the end of this section. 

To appreciate the following discussion, it must be recalled that the at- 
tenuation potential of an acoustic suppressor is dependent upon the modal cut- 
off ratio or angle of incidence of the sound entering the suppressor. Modes 
near cut-off (hitting the liner nearly normal) can be attenuated much more 
easily than modes of high cut-off ratio with mainly axial propagation in the 
duct. Roughly, the attenuation is inversely proportional to mode cut-off 
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ratio* Near cut-off modes radiate nearly normal to the engine axis and can 
therefore have a large Impact on fly-over noise, since the diatoms between 
the aircraft and the observer are minimized, Nearly axial propagating nolae 
usually travels a long distance and Is thus considerably attenuated when It 
reaches the observer, 

The angular location of the peak of the principal lobe of radiation for 
both Inlets (see reference 4, equation (33)) and exhaust ducts (equation (20) 
this paper) are shown as a function of mode cut-off ratio In figure 8. Notice 
that for an inlet at 0,4 Mach number the maximum angle of the principal lobe 
peak is about 66° at cut-off (6j) ** 1). All of the radiation at this angle and 
larger angles is from near cut-off modes, Thus the direct sideline radiation 
is quite easily attenuated. For an exhaust duct, however, the radiation at 
the typical radiation peak of 60° to 70° (see figures 5, 6 and 7) and forward 
towards the direct sideline Is composed of principal lobe radiation associated 
with quite large mode cut-off ratios. This sound will be attenuated much less 
by a suppressor. Thus, it appears that aft duct sound may be inherently more 
difficult to effectively suppress than that emanating from the inlet. Another 
consideration, on the favorable side, is that typically there Is more area and 
greater duct length- to-height ratio available in an exhaust duct than in an 
inlet. 

It is the author's opinion that there is an additional factor working to 
give less apparent suppression in the exhaust suppressor. Due to the abrupt 
nature of the exhaust nozzle and the near normal (to the duct axis) direction 
of the wave front vector for even moderate cut-off ratios, a strong reflection 
may occur for modes near cut-off. if this is true, then the duct modes which 
would be absorbed quite easily by the duct suppressor (near cut-off) never are 
observed in the far-field. It is thus possible that very small far-field at- 
tenuations could be observed even though large in-duct attenuations are ob- 
served in the exhaust duct. 


CONCLUDING REMARKS 

The far-field radiation emanating from the end of an exhaust duct has been 
studied using both approximate and exact methods. This initial study simulates 
the radiation of internally generated noise (as from a turbofan) from a circular 
exhaust duct with a resulting jet which differs in velocity from the surround- 
ing medium across a slip layer, The approximate method starts with the rela- 
tively simple single mode flanged duct radiation equation for zero flow. A se- 
ries of angle transformations adds flow and thus convective effects and also 
refractive effects across the slip layer. The exact method is an available 
single mode Wiener-Hopf radiation model for a circular duct. The approximate 
radiation pattern is a very accurate reproduction of the exact results includ- 
ing the radiation principal lobe shape and location as well as the location of 
the zeroes in the radiation pattern. An approximate multimodal radiation pat- 
tern was also derived which provides a simple expression for radiation with 
equal acoustic power per mode, This also agrees well with an exact equal 
acoustic power per mode radiation pattern generated using the Wiener-Hopf solu- 
tions. Experimental data of narrow band tone noise from static tests were com- 
pared to the multimodal radiation pattern. In most cases, for the blade pas- 
sage fundamental, the data agree qualitatively with the theory. For the second 
harmonic if is apparent that some dominant radiation is present from a limited 
number of modes. Finally, some considerable differences between inlet and ex- 
haust radiation became apparent while working with the approximate radiation 
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theory. It is pointed out that possibly the exhaust noise in the far-field 
is inherently more difficult to attenuate than the inlet noise using duct sup- 
pressors . 
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TABLE I. - DESIGN CHARACTERISTICS OF FANS A, B, AND C 


Characteristic 

Fan A 

Fan B 

Fan C 

Number of rotor blades 

40 

26 

26 

Number of bypass guide vanes 

90 

60 

60 

Rotor inlet tip radius, m 

0.9315 

0.9315 

0.8675 

Outside radius of nominal 
bypass nozzle,® m 

0.827 

0.827 

0.797 

Inside radius of nominal 
bypass nozzle,® m 

0.556 

0.556 

0.562 


a A£ter trimming. 


TABLE II. - AERO-ACOUSTIC CHARACTERISTICS OF FANS A, B, AND C 


Characeristlc 

Fan A 

Fan B 

Fan C 


Takeoff 

Approach 

Takeoff 

Approach 

Takeoff 

Approach 

Fan rotational speed, rpm 

3189 

2126 

3209 

2135 

4671 

3117 

Fundamental blade passage 

2126 

1417 

1390 

925 

2024 

1350 

frequency, Hz 







Frequency parameter,® q 

11.7 

7.77 

7.61 

5.06 

10.3 

6.88 

Bypass flow, kg/s 

379.2 

252.9 

381.3 

253.8 

372.7 

248.6 

Bypass exhaust Mach number a 

0.708 

0.472 

0.711 

0.473 

0.797 

0.532 




















































RELATIVE SOUND PRESSURE LEVEL, dB 




Figure 3. - Comparison of principal lobe peak and sldelobe zero locations of approxi- 
mate and exact solutions (% ■ 0.6( Mf ■ 0). 









